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The cosmic abundance of a long-lived charged particle such as a stau is tightly constrained by 
the catalyzed big bang nucleosynthesis. One of the ways to evade the constraints is to dilute those 
particles by a huge entropy production. We evaluate the dilution factor in a case that non-relativistic 
matter dominates the energy density of the universe and decays with large entropy production. 
We find that large Q balls can do the job, which is naturally produced in the gauge-mediated 
supersymmetry breaking scenario. 

I. INTRODUCTION 



Gauge-mediated supersymmetry (SUSY) breaking (GMSB) [l[ is appealing since the problems associated with the 
dangerous flavor-changing processes and CP violations are elegantly solved. In this scenario, the gravitino mass to 3 / 2 
is lighter than the weak scale and therefore the gravitino is most probably the lightest supersymmetric particle (LSP). 
■ If the i?-parity is conserved, the gravitino LSP is absolutely stable and can be a good candidate for dark matter 
' (DM) [3, d, 0,H, S3- In addition, such a scenario with a light gravitino may lead to spectacular collider signatures 
especially if the next-to-lightest superparticle (NLSP) is a electrically charged particle such as a stau. The stau 
NLSP has a quite long life time, since the decay rate is suppressed by the Planck scale Mp(~ 2.4 x 10 18 GeV). Such 
longevity may enable us to measure the Planck scale at collider experiments and therefore to test supergravity [8[, if 
P -1, the gravitino mass is O(10) GeV. 

However, the existence of such long-lived charged particles can jeopardize the success of the big bang nucleosynthesis 
(BBN). Since the stau decays during or after BBN, the energetic decay products may alter the primordial abundance 
of the light elements Furthermore, it was recently found that the negatively charged (long-lived) particle can 

form bound states with nuclei (e.g., 4 He— f ~) , which catalyze the nuclear reactions and substantially change the 
abundance of the light elements such as 6 Li 11]. The detailed discussion based on the catalyzed BBN (CBBN) limits 
the life time of the stau as < 10 3 sec, assuming the thermal relic abundance for the stau. Taking the upper bound at 
face value, we are led to smaller m 3 / 2 and/or heavier nif, which are challenging from the experimental point of view 
\^ or on the basis of naturalness. 

^vq ■ To alleviate the BBN constraints, several solutions (e.g., i?-parity violation |l2j) has been proposed. Among them, 
it is simplest to assume that there is late-time entropy production to dilute the stau abundance. According to Ref. [Hj], 
the necessary dilution factor is A = (300 — 600) x (wf/lOOGeV) for m 3 / 2 = 10 GeV, where m? denotes the stau mass. 
However, successful late-time entropy production is not so easily achieved as one might think of. The reason is as 
follows 14]. To be concrete, we assume that the gravitino is LSP while the stau is NLSP, and that a scalar field X with 
even i?-parity produces large entropy. Then the scalar mass mx must be lighter than m?, since the supersymmetric 
partners of the standard-model (SM) particles, if kinematically allowed, are generically produced from the decay of 
X [Tolfilil]. In addition, since mx is lighter than the stau mass, X must be more strongly coupled to the SM particles 
than the gravitational interactions in order to decay before BBN. On the other hand, the fermionic partner of X 
must be heavier than nif, because we have assumed that the stau is NLSP. Therefore the mass spectrum must satisfy 
mx < rrif < niy, where m x is the mass of the fermionic partner of X. To realize such spectrum, however, one needs a 
partial cancellation between the SUSY mass and the soft SUSY breaking mass. Thus, it is difficult, if not impossible, 
to naturally induce a late-time entropy production in the set-up with the gravitino LSP and the stau NLSP. 

In this article, we show that the successful entropy production can be achieved by the decay of the Q ball in the 
minimal supersymmetric SM (MSSM). In our scenario, a large charge Q naturally makes the effective mass smaller 
than m? and it also explains the longevity of the Q ball. Since the fermionic partners are nothing but the SM particles, 
there is no constraint on the fermion mass. In addition, a right amount of the baryon asymmetry can be generated 
by the Affleck-Dine (AD) mechanism [13, El- 

In the next section, we provide a brief review on the Q balls in GMSB. We estimate the dilution factor in both cases 
when the Q ball dominates the universe after and before the freeze-out of the stau in Sec. Mil In Sec. IIV1 we show 
that the Q-ball decay can dilute the stau by the desired amount. In Sec. [V] we give our conclusion and discussions. 



' This ameliorates a possible gravitino production from the scalar field X |16|1 . 
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II. Q BALL IN GMSB 



A Q ball is a non-topological soliton of a complex scalar field $ given by the minimum energy configuration with 
a fixed U(l) charge Q |l9|. One of the conditions for the Q balls to be formed is that the scalar potential with a 
U(l) symmetry is shallower than the quadratic potential at large field value. The Q balls are known to be formed 
associated with the scalar dynamics of the MSSM fields, especially in connection with the AD mechanism [UUli]. In 
addition, since the Q balls generically have a very long life time, they can play important roles in cosmology. 

In MSSM, there are many flat directions composed of some combination of squarks, sleptons and Higgs bosons. 
Alon g th e flat directions, both .F-term and D-term potentials vanish in the exact SUSY limit at renormalizable 
level [l8l [20j . They are lifted by the soft SUSY breaking effects, non-renormalizable operators, and finite temperature 
effects. In the AD mechanism, one of the flat directions (denoted by $) is assumed to have a large field value during 
inflation *. After inflation, $ starts to oscillate when the Hubble parameter becomes comparable to the mass of At 
the same time, $ acquires the baryon (and/or lepton) asymmetry, due to non-renormalizable baryon-(lepton-)number 
violating operators that are effective only at large field values. The scalar potential of $ has an approximate U(l) 
symmetry corresponding to the baryon and lepton symmetries that are conserved at low energy effective theory, i.e., 
MSSM. As we will see below, since the scalar potential is shallower than the quadratic potential above the messenger 
scale in GMSB, $ experiences spatial instabilities and deforms into Q balls, where the charge Q corresponds to the 
baryon and/or lepton numbers. 

As mentioned above, the scalar potential is lifted by the SUSY breaking effects, non-renormalizable operators, and 
finite temperature effects. In GMSB, the scalar potential above the messenger scale is given by 

/ ld>l 2 \ 2 / I5>l 2 \ A 2 |ff>| 2 (™- 1 ) |<j>| 2 

U($) ~ M% (log J + c g ml /2 \l + K log LL j |$| 2 + AJ__ + CtT 4 lQg IJ_ _ CHH 2^ ?y (1) 

where m 3 / 2 is the gravitino mass, Mp ~ 2.4 x 10 18 GeV is the reduced Planck scale, and we omit the baryon- (lepton- 
)number violating operators here. The first term comes from the GMSB effect above the messenger scale Ms [HI]. 
Mp and Ms are related to the F and A components of a gauge-singlet chiral multiplet S in the messenger sector as 

M ' F= (4^ K2(F5>2 ' M * = ^), (2) 

respectively, where g collectively stands for the SM gauge couplings, and k denotes the Yukawa coupling constant 
between S and the messenger fields. In general, Mp could be in the range 10 3 GeV < Mp < O.I^/m 3 / 2 Afp ~ 5x 1 8 
GeV for to 3 / 2 = 10 GeV. The second term of Eq. fl} comes from the gravity- mediated SUSY breaking, and the 
coefficient c g is of the order unity. Here the one-loop correction is included, and K is negative with \K\ = 0.1 — 0.01 
for most of the flat directions. Since the gravitino mass is relatively small in the gauge mediation, this term is 
effective for a large field value of $. The third term in the potential is due to a non-renormalizable interaction in the 
superpotential of the form Wnr = A$ n /(nAC~ 3 ) with n > 3, where A is a coupling constant. The fourth term is a 
two- loop thermal correction to the potential [22j . The coefficient cp can be both positive and negative depending on 
which flat direction we choose [23[, and the absolute value is roughly given by fp — O(0. 1), where we define \cp\ = f T - 
The last term is a Hubble-induced mass term, which stems from the quartic coupling in the Kahler potential between 
the flat direction and the inflaton [l8| . This term is absent after the reheating of the inflaton. Because of this term, 
the potential has a minimum at a large field amplitude during inflation, and the $ field is trapped there and it serves 
as the initial condition for the later dynamics To be concrete, we assume throughout this paper that the minimum 
is given by the balance between the Hubble-induced mass term and the non-renormalizable operator: 

i 

11 '""m p , (3) 



XM P 



where 4> = v2|$|< The flat direction $ traces the minimum until it starts to oscillate. 

After inflation, the flat direction starts rotating, experiences spatial instability, and deforms into Q balls [25|, [26|, 
US HH Hi| ■ The properties of the Q ball are well known. The charge Q is determined by the amplitude of the $ field 



t Actually, multiple flat directions can have large expectation values simultaneously, but we do not consider this possibility for simplicity. 
§ The same effect may be realized by large enough Hubble-induced A terms [24ll , 
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at the onset of the oscillations, 4> OS c- If the potential at cj) = <fi osc is dominated by the first or fourth term in Eq.([T]), 
the charge of the Q ball, which is called the gauge-mediation type, is determined as [29] 



where M(T) is defined as 



13 (#)) ' (4) 



M(T) = { f r %> If™ (5) 

1 /t^ osc for f T T osc > M F w 



The subscript "osc" denotes that the variable should be evaluated at the onset of the oscillation of <E>. The numerical 
coefficient (3 is given by [U 

J 6 x l(T 4 e for e > 0.1 , , 

" " [6x 10~ 5 for e < 0.1 ' W 

where e (< 1) denotes the ratio of the baryon number density to the number density of The size and mass of the 
Q ball, and the effective mass of the field inside the Q ball (i.e., the mass per unit charge) are written respectively 

as m n 

I 1 4a/27T S i— 1 

7Im(T) QI ' Mq*^-M{T)Q*, co Q ^V27rM(T)Q-. (7) 

On the other hand, if the potential at <f> = 4> osc is dominated by the second term in Eq.([T]), the gravity- mediation 
type of Q ball is formed, whose charge is given by [1^, [3(| 

Q = {?(^)\ (8) 



where 



V m 3/2 



R , . 6 X 10~ 3 e for e > 0.1 , , 

P ~ i 6 X 10- 4 for e < 0.1 ' 1 j 



The size and mass of the Q ball, and the effective mass of the field inside the Q ball (i.e., the mass per unit charge) 
are written respectively as |26l l29l |30T ] 

R Q~T77a ' M Q- m 3/2Q, UQ-m 3/2 - (10) 

\K\2m 3/2 

The Q ball can decay if the mass per unit charge loq is larger than the decay products that carry the same charge. 
For example, if Q is the baryon number, the lightest particle with baryonic charge is a nucleon whose mass is ~ 1 GeV. 
Therefore, such Q balls with uiq > 1 GeV can decay into the nucleons (perhaps together with other lighter particles 
such as 7r- mesons). Since the decay can proceed only from the surface of the Q ball, the decay rate is bounded from 
above. For the MSSM Q ball, the rate is saturated and given by [3l| 

-=Q 4 for gauge-mediation type 

r — :Q for gravity- mediation type 

24n\K\ 

Therefore, the decay temperature of the Q ball is calculated as 

_ i 



10 MeV g*, 4 I 1Q7 j I J for gauge-mediation type 

KY\~* ( 7123/2 Y (-9. 



5 MeV ( — J {^q^ ) I 10 s ' f ° r g ravit y- mediation ^P 



(12) 



where g* q counts the relativistic degrees of freedom at the Q-ball decay, and we define g* q = c?*^/ 10.75 in the second 
equality. 
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FIG. 1: Evolution of the energy densities of the radiation (dashed (red)), the Q-ball (solid (blue)), and the radiation from the 
Q-ball decay (dashed dotted (green)). 



Let us now estimate how much the abundance of the stau is diluted in a situation that the Q ball releases entropy 
after it dominates the energy density of the universe. We define the dilution factor A as 



where the left-hand side represents the ratio of the stau number density n? to the entropy density s in the presence of 
the entropy production due to the Q-ball decay, while the abundance of the stau on the right-hand side is estimated 
without the entropy production. Note that, although we consider the Q-ball decay, our arguments in this section 
can be applied to any scenario that non-relativisitc matter dominates the universe and decays with large entropy 
production. 

The stau abundance can be diluted if the Q-ball decay takes place after the freeze-out of the stau. The dilution 
factor depends on the thermal history. We assume that the universe is radiation-dominated before the Q balls start 
to dominate the energy density of the universe. Let us define T eq , Tf a , and To as the temperatures when the Q-ball 
energy density becomes equal to the radiation density, the stau freezes out, and the Q ball decays, respectively. (We 
will use such notation that the subscripts eq, fo, and D denote that the variables should be estimated at T — T eq , Tt , 
and Td, respectively.) In Fig. [T] we sketch the evolution of the energy densities of the radiation, the Q ball, and the 
radiation produced by the Q-ball decay. 

We now consider the following cases: (i) T eq < Tf and (ii) T eq > X/ . In the case (i), the situation is very simple, 
since the Q balls change the evolution of the universe after the freeze-out of the stau. We thus obtain 



III. ESTIMATE OF DILUTION FACTOR 




(13) 




(14) 



The dilution factor is calculated as 




T 



(15) 
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where pq is the energy density of the Q balls. 

In the latter case (ii), the stau decouples from thermal equilibrium when the Q ball is dominating the universe. 
Therefore one needs to know the dependence of the stau number density on the Hubble parameter at the freeze-out. 
Since the freeze-out takes place when the annihilation rate becomes comparable to the expansion rate, we obtain 

nrjo ~ fe, (16) 
(crv) 

where (crv) is the thermally averaged cross section for the stau annihilation processes. The freeze-out temperature 
becomes larger than that in the usual radiation-dominated universe, since the energy density at the freeze-out is 
higher in the presence of the Q balls. However the change in Tf D is not significant, because the stau decouples when 
it is non-relativistic, so the freeze-out temperature is rather insensitive to the change in the cosmic expansion rate. 
Thus, we simply assume Tf ~ m^/20 in the following Such an approximation is not essential to our arguments. 
Then we can estimate A as 



") (-) (-) 



PQ) D \ n f/ fo 
Tfo I Hf 



thermal 



(17) 



where we use pq ~ 3Hf a Mp at the freeze-out, and H^* denotes the Hubble parameter at the freeze-out in the 
absence of the Q-balls. 

The relation betwen the Hubble parameter, Hf a , and the freeze-out temperature, Tf Q , depends on whether the 
newly created radiation from the Q-ball decay dominates over the radiation that was present from the beginning, at 
the freeze-out of the stau. Let Tt r denote the transition temperature at which both radiation components become 
comparable to each other. (See Fig. Q]) In the case of T/ D < T tr , we have 

where we use T tr ~ (TeqTfS)i and H fo ~ T^T^Mp 1 . On the other hand, in the case of Tf > T tr , we have 
H fo ~ (Tf T eq )i /M P . Substituting this into Eq. (H)), we obtain 

A ~ gV^kll . (19) 
Td 



In summary, we find that the dilution factor is given as follows depending on when the freeze-out takes place: 

(Case A : T eq < T fo ) 

(Case B : T tr < T fo < T eq ) . (20) 
(Case C : T D < T fo < T tr ) 



T 

1 eq 

Td~ 

(TfoT eq )i 
T D 
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Tfo 
T D 

Note that the decay temperature should lie in the range of 5 MeV < Td < Tf Q [32j and it must satisfy Td < T eq , in 
order to dilute the stau by the entropy production. 

In the next section, we will determine the Q-ball-radiation equality temperature T eq by considering the formation 
and the dynamics of Q balls, in order to evaluate the dilution factor A. 



■ To include the change in Tf a , let us define 7 = Tf /T'j o , where T^ o denotes the freeze-out temperature for the modified cosmic 
expansion. It is estimated as ■yg ~ 1 — (T/2m,f ) log(T eq /Tf ) for the case B, and 70 ~ 1 — (2T/m,f ) log(Tf /Trj) for the case C (See 
Eq J20l l and Fig. [TJ. The freeze-out temperature does not change in the case A. As long as we stick to A ~ 10 3 , the minimum values of 
IB, C & re given by 73 ~ ■yc ~ 0.7, and our approximation seems to be valid. To take account of this change in Tf a , one has to multiply 
A in Eq. J20} by 7" 1 in the cases B and C with Tf a replaced with Tf - Then, one can see that A becomes larger by a factor of ~ 4 at 
most. Note that including the effect always increases the dilution factor. 
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IV. ENTROPY PRODUCTION BY THE Q-BALL DECAY 

The Q-ball formation and the subsequent thermal history depend on the scalar potential at the onset of the 
oscillations. Since the large charge Q ball is necessary for a long lifetime, the field amplitude at the onset of the 
oscillation should be very large. In that case, the zero-temperature potential would be dominated by the gravity- 
mediation term (the second term in Eq.fl])). As we will see below, the thermal corrections are negligible for the 
parameters we adopt in the following analysis. 

The flat direction starts to oscillate when 3H ~ m 3/ / 2 . Here we simply assume that it takes place when the inflaton 
oscillation dominates the energy density of the universe. Then the ratio of the energy densities of the Q ball and the 
inflaton at the reheating is given by 

PQ 



Pvnf 



RH P m f 



2 m 3/2^oac 3 / <?W X " 



3(m 3/2 /3) 2 M 2 2 V M 



After reheating the ratio evolves as cx T 1 , so the Q-ball-radiation equality temperature is obtained as 



O / ± \ 2 / rji \ / I \ 2 

o I (pose \ - „ , r / J RH 



(21) 



Te - r -2^J ~ 5 GeV U x 10 7 GeV J U x 10 14 GeV ) ' (22) 
Meanwhile the decay temperature is calculated from Eq. (fl2|) as 

Td -{-9Q-) \-24^K\) P ^;- 5MeV Uxl0 14 GeV) llOGevJ ' (23) 

where \K\ = 0.01 and g* q — 10.75 are used. Since the freeeze-out temperature is Tf a ~ 5 GeV for rrif = 100 GeV, 
the dilution factor is estimated as in the case A: 



A ~^i~l0 3 f ^ ) ^ ) ( W3 / 2 1 ' 

T D I 8 x 10 7 GeV/l 5 x 10 14 GeV/ 1 10 GeV ; ' 



In this case we have Q ~ 10 24 . Notice that Tv shown in (|2"3")l can marginally satisfy the BBN constraints [32]. For 
slightly larger TO3/2 or smaller 4> osc , one can have a large enough Tp that safely satisfies the BBN bound, keeping 
the dilution factor A large enough. Since the initial amplitude of the flat direction is determined as in Eq.([3]), it will 
be realized for the n = 6 direction (LLe or udd) with A ~ 0.006, or the n = 7 direction (dddLL) with A ~ 30. The 
coefficient of the thermal logarithmic corrections will be negative for these directions. If the thermal correction to 
the potential dominates over the gravity-mediation term, it will spoil the above scenario because the $ field will be 
trapped by the negative thermal logarithmic potential, and so, it cannot be released for a long time. In order to avoid 
such a situation, we must impose a condition that the thermal logarithmic correction is negligible at the onset of the 
oscillations: j^tsc < \ m \ii^ 2 osc **■ Using T 4 SC ~ 0.5 H osc M p Tr H , it is rewritten as 



> ^ (*L) ^ 2 x 10" GeV (£) 2 (- T ™\ k , (25) 



s/l \m 3/2 J \0.lj \8 x 10 7 GeV/ VlOGeV 



which is satisfied in the above analysis. 

Finally, we comment on the last moment of the Q-ball decay. As the charge becomes small, the gravity-mediation 
type Q ball gradually deforms into the gauge-mediation type one. Therefore, at a certain point, the mass per unit 
charge of the Q ball may exceed the stau mass, ojq > m,f, which implies that the stau can be produced. In order to 
suppress the stau abundance as Yf < 10~ 17 , we must impose Mf ^$ 10 5 GeV. This is because Yf ~ -Bf(Tb/m3/2) 
and the branching ratio is estimated as Bf ~ Q C r/Q where Q cr ~ (Af^/mf) 4 . Such a small value of Mp is realized 
in a model where the Yukawa coupling k is suppressed as in the case of the composite S field. 



The thermal logarithmic term may not appear in such situation that the dilute plasma before reheating is suppressed because of e.g.. 
the small mass of the inflaton and/or the existence of appropriate multiple fiat directions. 
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V. CONCLUSION 



We have shown that the Q-ball decay can produce large enough entropy to dilute the cosmic abundance of a long- 
lived charged particle such as a stau. Since the Q balls are composed of the MSSM particles, our scenario is minimal in 
some sense. Successful late-time entropy production is not so easily achieved as one might think of: the particle that 
produces entropy should have a smaller mass than the stau and a long lifetime, and no unwanted particle production 
should occur. In our scenario, a large charge Q naturally makes the effective mass smaller than the stau mass, and 
the decay products are just the SM particles. It also explains the longevity of the Q ball. The large Q ball can be 
naturally produced in the dynamics of the flat direction. 

In addition, we have derived analytically the dilution factor for the cases that the stau freezes out both before and 
after the Q ball starts to dominate the universe. 

Lastly let us briefly discuss the baryon asymmetry and the dark matter. Due to the late-time entropy production, 
the baryon asymmetry and the gravitino dark matter are also diluted. One way out is to over-produce both by the 
amount of the dilution beforehand. Another is to create them after the entropy production. As for the gravitino 
dark matter, the thermal production does not suffice since the reheating temperature cannot be too high due to 
the constraint ([25]) . Therefore, one has to rely on the non-thermal production @, l33| . On the other hand, one 
can obtain a right magnitude of the baryon asymmetry from the Q-ball decay itself, making use of the Affleck-Dine 
baryogensis. If the flat direction has the baryon number, it must start the oscillations with suppressed angular motion 
(i.e., e ~ O(10~ 7 )). On the other hand, for the leptonic direction, the lepton charges evaporated from the Q ball 
before the electroweak phase transition are partially converted into the baryon asymmetry through the sphaleron 
processes. In this case, e ~ O(10 -3 ) is necessary to have a right abundance of the baryon asymmetry. 
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